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Abstract 

I consider a theoretical description of recent experiments on doping the spin-Peierls 
compound CuGeC>3 and the Haldane gap compounds PbN^X^Os and Y2BaNiOs. The 
effective theory is the one of randomly distributed spin-1/2 moments interacting with 
an exchange decaying exponentially with distance. The model has two phases in the 
(doping, interchain coupling) plane: (i) A Neel ordered phase at small doping; (ii) A 
quantum disordered phase at larger doping and small interchain interactions. The spin- 
Peierls compound CuGeC>3 and the Haldane gap Nickel oxides PbNi2V20s and Y2BaNiOs 
fit well into this phase diagram. At small temperature, the Neel phase is found to be 
reentrant into the quantum disordered region. The Neel transition relevant for CuGe03 
and PbNi2V208 can be described in terms of a classical disordered model. A simplified 
version of this model is introduced, and is solved on a hierarchical lattice structure, 
which allows to discuss the renormalization group flow of the model. It is found that the 
system looks non disordered at large scale, which is not against available susceptibility 
experiments. In the quantum disordered regime relevant for Y2BaNiOs, the two spin 
model and the cluster RG in the ID regime show a power law susceptibility, in agreement 
with recent experiments on Y2BaNiOs. It is found that there is a succession of two distinct 
quantum disordered phases as the temperature is decreased. The classical disordered 
model of the doped spin-1 chain contains already a physics relevant to the quantum 
disordered phase. 
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1 Introduction 



Doping quasi one dimensional (ID) antiferromagnets with a spin gap has become experimentally 
possible since the discovery of several inorganic quasi ID oxides. One of these compounds is CuGe03 
having a spin-Peierls transition at T$p — 14 K [||. Below T$p, the spin-phonon coupling induces a 
dimerization of the lattice, and the opening of a gap in the spin excitation spectrum. The Haldane gap 
in spin-1 chains in another example of a spin gap state in low dimensional magnets @. Two inorganic 
spin-1 Haldane gap antiferromagnets have been discovered in the recent years: (i) PbN^A^Og having 
a spin gap ~ 28 K |J. (ii) Y 2 BaNi05 having a spin gap ~ 100 K M. The spin-Peierls compound 
CuGe03 and the two Nickel oxides PbNi2V20s and Y2BaNiOs can be doped in a very controlled 
fashion. Substituting the magnetic Cu sites (having S = 1/2) of the spin-Peierls compound CuGeOs 
with a variety of ions (Ni || - a spin-1 ion -, Co Q - a spin-3/2 ion -, Zn |7|, ||, [9], [II], 11] or 



Mg [12] - non magnetic ions -), or substituting the Ge sites with Si [T^] leads to the formation of 
an antiferromagnetic phase (AF) at low temperature. Moreover, in CuGe03, there is AF long range 
order even with an extremely weak concentration of Zn impurities [14]. On the other hand, the 
two Nickel oxides PbNi2V208 and Y2BaNiOs have been the subject of an important experimental 
interest recently. It has been shown that substituting the spin-1 Ni sites of the PbNi2V20s compound 
with Mg - a spin-0 ion - leads to AF long range order. In the Y2BaNiOs compound, the Ni sites 
can be substituted with Zn or Mg - non magnetic ions -. In this case, no sign of AF long range 
order has been reported, even at extremely low temperature [15, 17]. Instead, it has been found 
that the susceptibility has a power-law temperature dependence [18|. Experiments therefore show 
that doping quasi ID antiferromagnets with a spin gap can lead to very different situations: either 
antiferromagnetism, or a power-law susceptibility without AF ordering. The purpose of the present 
article is to describe these experimental observations in a unified theoretical framework, and provide 
a detailed theoretical analysis of the different phases of the model. 

On the theoretical side, a lot of efforts have been devoted to understand the behavior of random 
spin chains. The theoretical tool usually used to study these models is the cluster renormalization 
group (RG) |ll|, which is a perturbation theory in the inverse of the strength of the strongest 
exchange in the chain, and leads to a certain number of exact results at low temperature because 
the exchange distribution becomes extremely broad. A lot of different models have been solved using 



this approach. For instance: the Ising chain in a transverse magnetic field [20], the random spin- 



1/2 chain |2l[], the spin-1/2 chain with random ferromagnetic and antiferromagnetic bonds [22], the 
dimerized chain with random bonds [23], the disordered Haldane gap chain [24, 25 1. These studies 
have revealed that ID disordered magnets can be controlled by several types of Griffiths phases: (i) 
the random singlet phase with a diverging susceptibility and algebraic correlation; (ii) the "weakly 
disordered" phase with a diverging susceptibility and short range correlations. 

An important question is to understand the relation between the available experiments on quasi 
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one dimensional oxides and the available theories of disordered ID magnets. This type of approach 
followed in the present article incorporates realistic constraints such as interchain couplings and a 
finite temperature. 



The starting point of such our description has been already established in previous works []2q, 27 



28, 29, 30 1 . In the doped spin-Peierls systems, non magnetic impurities generate solitonic spin-1/2 



degrees of freedom distributed at random with a concentration x. The solitons are confined close to 



the impurities because of interchain interactions [30, 31, 33, 34, 35, gg, 37], and interact with the 
Hamiltonian 

H = ^2JijSi.Sj. (1) 

(id) 

The exchange between two spin-1/2 moments at positions (xi,yi) and (xj,yj) is mediated by virtual 
excitations of the gaped medium and therefore decays exponentially with distance: 

= I ) " • " " Aexp ^(^-^) 2 +(^-^) 2 j , (2) 
where £ x and £ y are the correlation lengths in the direction of the chains and perpendicular to the 



chains respectively [28, 29. |30| ]. The form Eq. Q of the exchange incorporates a correlation length 
in the transverse direction shorter than in the longitudinal direction (£y = £ z /10, and £ x . ~ 10 
in CuGe03 |39|]). The exchange Eq. || is staggered because the dimerized pattern propagates 
staggered antiferromagnetic correlations. The Hamiltonian Eqs. [l], |2] is therefore strongly disordered 
but unfrustrated. 

Now, the model relevant to describe doping in a Haldane gap system is almost identical. It is 
well known that an impurity in a Haldane gap chain generates two "edge" spin-1/2 moments: one 
at the right and one at the left of the impurity site [40, 41, |42]| . The spin-1 chain can be thought in 



terms of a Valence Bond Solid (VBS) |43|]. Introducing a paramagnetic site breaks two VBS bonds, 
therefore resulting in two "edge" spin-1/2 moments. At energies far below the Haldane gap, only 
these edge moments are the relevant degrees of freedom (see Fig. [I]). The two edge moments in 
the same unit (i.e. generated by the same impurity) interact with a ferromagnetic exchange —J2 
originating from the coupling to neighboring chains, with therefore the same order of magnitude as 



the interchain interaction: J2 ~ Jj_ [44]. The edge moments belonging to different units are coupled 
by the staggered exchange Eq. ||. 

We find that, depending on the doping concentration and interchain interactions, the model 
Eqs. [l], H] has two regimes: a Neel ordered region and a quantum disordered region. In the Neel 
ordered region, relevant for CuGeOs and PbN^A^Os, quantum mechanics plays little role, and we 
are lead to replace the spin variables in Eq. |l], |2| by classical Ising spins. We propose here that this 
type of disordered Ising model is equivalent to another type of disordered Ising model, and solve 
the hierarchical lattice version of the latter. Using this treatment, we can compute the renormalized 
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J 2 (-)'a exp(-l/ j 2 



Figure 1: The low energy effective model of the doped spin-1 system. A paramagnetic impurity 
(square symbols) generates a unit of two edge moments. The edge moments originating from the same 
impurity are coupled ferromagnetically by an exchange — J2 having the same order of magnitude as the 
interchain coupling. Two edge moment at a distance I are coupled by the exchange (— )'Aexp (— £/£)■ 



exchange distribution. In spite of a strongly disordered initial Hamiltonian (see Eqs. |l], ||), we 
find that at large scale, the problem behaves as if it were non disordered. This appears to be 
consistent with susceptibility experiments showing a well-defined transition even with a small doping 



concentration p4 |. 

In the quantum disordered region of the phase diagram, the physics is dominated by the formation 
of randomly distributed singlets. We show that the susceptibility has a power-law behavior, which 



turns out to be in agreement with existing experiments on Y2BaNiOs [18|. We also find the existence 
of two distinct "quantum disordered" phases. The high temperature "quantum disordered" phase 



appears to have been observed experimentally in Y2BaNiOs [15]. There appears to be another low 
temperature "quantum disordered" phase in which the edge spins in the same unit (see Fig. |l]) are 
frozen into spin-1 objects. Finally, we show that even in the quantum disordered regime of the model, 
part of the quantum disordered behavior is already contained in the classical disordered magnet. 

The article is organized as follows: the phases of the model Eq. |2| are discussed in section |2[ 
We show the existence of two phases: a Neel ordered region and a quantum disordered region. The 
nature of these two phases is next discussed in details in sections || and Concluding remarks are 
given in section || 



2 Phases of the model 

In this section, we use a phenomenological approach to derive the phase diagram of the model. 
The calculation of the relevant energy scales in the problem is based on the analysis of a two-spin 
model. There is a first temperature scale (being a fraction of the spin gap A) below which magnetic 
correlations start to develop inside the chains. There is a second energy scale T typ , equal to the 
typical exchange, associated to singlet formation. There is a third energy scale Tstoner associated to 
long range AF ordering. The behavior of the model depends strongly on whether Tstoner is larger or 
smaller than T typ . 
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Figure 2: Temperature dependence of the internal energy of the classical and quantum two-spin 
model, with A = 44.7 K, £ = 10, x = 0.01, relevant for CuGeOs. The temperature and energy are 
calculated in Kelvin. In the two-spin classical model, we considered two spins coupled by an Ising 
term TC = (J/2)criO"2. 

2.1 Onset of magnetic correlations 

Magnetic correlations start to appear inside the chains when the temperature is a fraction of the spin 
gap A. To show this, let us consider a simple model in which two spins at a distance I are coupled 
antiferromagnetically: Ti = J(Z)Si.S2- We use an exchange decaying exponentially with distance (see 
Eq. pn: J (I) = Aexp (—£/£), and a Poisson bond length distribution V(l) = rrexp (— xl). Rigorously, 
the spacing I is a discrete quantity, distributed according to a geometrical distribution. However, the 
physics will turn out to be controlled by the large-Z behavior and it is legitimate to consider I as a 
continuous variable, and replace the geometrical distribution by the Poisson distribution. 
The internal energy of the two-spin model reads 



Magnetic correlations start to appear in the two-spin model in the low temperature regime in which 





This expression can be expanded in T: 




(4) 



U(T) is linear in T. This regime appears when T is a fraction of A (see Fig. |2|). 
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2.2 Singlet formation 

To discuss in what temperature range is the physics controlled by the quantum mechanical ground 
state (being a singlet), we need to calculate the probability V S (T) to find the two spins in a singlet 
state at a finite temperature T = 1/(3. We have 

V(T) - fdlV(l) exp [3/3J(Q/4] 

T\* ^ A „,_ 1+x , exp(-u/4) 



= l-3sf 7 / u~ 1+ ^ ; ^ V —du, (6) 

7o exp(3u/4) +3exp(-u/4) v; 

where we used the dimensionless parameter u = j3J. The integral in Eq. ^ is dominated by the small 
exchanges and we have V S {T) ~ 1 - k(T/A) x ^, with numerical factor. We are lead to conclude 
that the ground state occupancy is close to unity below the energy scale Tt yp ~ A exp (— l/(x£)). Tt yp 
is nothing but the typical exchange, already identified in a previous work on the ID model [27J. 

The calculation of the ground state occupancy can be made even simpler by noticing that only 
the disorder configurations in which the exchange is larger than ~ T are in a singlet configuration. 
This leads to 



V S (T)~J xexp(-xl)dl ~ 1 - \ 



It is remarkable that the energy scale Tt yp arising from the two-spin model is identical to the one 



obtained previously from the exact solution of the ID effective model of the spin-Peierls chain [27|. 
This shows that a model with only two spins contains already the relevant physics. 

2.3 Neel ordered versus quantum disordered behavior 

At low temperature, correlations between the chains induce a long range ordering of the spin system. 
The simplest phenomenological description of long range ordering is provided by a Stoner model, 
already considered in Ref. p8| . In CuGeC>3, there is a succession of three regimes (see Fig. ||-(a)): 
(i) a paramagnet at high temperature, (ii) intrachain correlations develop when the temperature is a 
fraction of A (iii) interchain correlations give rise to long range antiferromagnetism below Tstoner = 
J±x£. 

In PbNi2V20§, the relevant parameters are A ~ 30 K, x ~ 0.02, J± ~ 1.1 K. We approximate 
the correlation length in the Haldane gap phase to be £ ~ 6. The true correlation length is expected 
to be larger than this value because PbN^X^Os is close to a transition to an Ising ordered antifer- 
romagnet H. We find Tg t0 ner = 0.13 K, T typ = 6 mK, showing that the same succession of regimes 
occur in PbN^A^Os and CuGeC>3 (see Fig. Q-(a)). This is compatible with existing experiments in 
PbNi 2 V 2 8 §. 

Therefore, in CuGe03 and PbNi2V20s, one has Tgtoner S> Tt yp . This implies that singlet forma- 
tion plays little role in the physics of the antiferromagnetic transition. The quantum two-spin model 
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Figure 3: In the models with 3D antiferromagnetism (a) T typ is far below Ts toner = J±x£, with 
therefore a classical transition to an AF phase. This is the case for PbNi2V208 as well as CuGeC>3. 
In the models with a quantum disordered ground state (b), Jj_ is far below T typ , which occurs in 
Y2BaNi05. In this case, the classical paramagnet has a cross-over to a ID quantum disordered 
state with formation of random singlets. Below J±, singlets are formed among spins belonging to 
different chains. Below Tg to , there is a transition to a reentrant antiferromagnetic phase. Intrachain 
correlations start to play a role when the temperature is below a fraction of the spin gap A. 



can then be well mimicked by the classical two-spin model. For instance, the internal energy of the 
classical and quantum two-spin models have an identical temperature dependence (see Fig. |2|). This 
indicates that one might expect to obtain a reasonable description of antiferromagnetism in CuGeOs 
and PbNi2V20g on the basis of a classical model, which we analyze in details in section ||. 

Now, the situation is different in Y2BaNiOs, where one has A ~ 100 K, £ ~ 6, J± ~ 0.3 K, and 
x ~ 0.04 Q. We find T typ = 1.6 K, and T St oner = 0.07 K. What is new compared to CuGe0 3 and 
PbNi2V20g is that a well defined quantum disordered regime is present below Ttyp- I n between J± 
and Ttyp, there is singlet formation in the chain direction, and below Jj_, the singlets develop in the 
transverse direction. Below the energy scale Tt yp , the staggered susceptibility is well described by 
the one of the random spin-1 chain [E4L E5fl 



X {T) ~ xT a ~ L /T° 



typ' 



(7) 



Using a cluster RG calculation, we will calculate the susceptibility in section 4.2 and show that it has 



indeed a power-law temperature dependence. The Stoner criterion leads to the ordering temperature 



rp/ 
Stoner 

T typ 



(J±x£X 

V T typ J 



This shows that the quantum disordered model transits to a reentrant antiferromagnetic ground state 
below T'g tonor . The different regimes of the model are shown on Fig. §-(b). 
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Figure 4: Phase diagram in the plane Jj_/A). The solid line is Jj_/A = l/(x£)exp [— l/(x£)], 
delimitating the transition from the Neel ordered region to the quantum disordered region. Various 
experimental systems have been reported on the diagram: +, x: Y2BaNiOs [18, 17] (being quantum 
disordered); *: CuGe03 |L4|] (being antiferromagnetic) ; □: PbN^X^Os |§ (being antiferromagnetic) . 
The values of J± and £ have been taken from these references. This phase diagram is valid above 
the temperature scale 2g toner . Below T stoner , the antiferromagnet is reentrant inside the quantum 
disordered phase. 



The existence of two classes of models is best summarized by calculating the ratio 

- c 1 



Tstoner ^-1- _/■,!/ (xf) 



^typ A 

which can be smaller or larger than unity, controlling whether the model has a Neel transition at 
Tstoner or is in a quantum disordered regime. These two behaviors, as well as a comparison between 
the model and existing experiments on CuGe03, PbNi2V20s and Y2BaNiOs, have been reported on 
the phase diagram on Fig. ||. 

3 Nature of the antiferromagnetic transition 
3.1 Motivation of the hierarchical lattice study 

In low doping experiments in CuGeC>3, Manabe et ah have measured the doping dependence of 
the Neel temperature, and found that the experimental data in the range x > 0.1% could be well 
fitted by the behavior T\r ~ Aexp (—B/x) J14[|. This suggests that there is no critical concentration 
associated to antiferromagnetism. As already pointed out in Refs. [28, 29], the model Eqs. |], [2| 



shows sings of compatibility with these experiments. The main unsolved question raised by the 
experiments by Manabe et al. jL4|] is to determine whether the maximum in the susceptibility is 
really the signature of an antiferromagnetic phase transition with a diverging staggered correlation 
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(a) (b) 

Figure 5: The hierarchical diamond used in the Migdal-Kadanoff renormalization. (a) A hierarchical 
lattice with 3 generations, (b) A hierarchical lattice with 2 generations. One RG transformation 
consists in decimating the sites at the deepest into the lattice, i. e. for instance transforming the 
lattice (a) into the lattice (b). 

length. In fact, the susceptibility experiments give no information about the existence / absence 
of a diverging correlation length, and there are are no neutron experiments with a concentration 
of impurities of order ~ 0.1%. On the theoretical side, we used previously several approaches to 
describe the nature of the antiferromagnetic phase of the model Eqs. |l], ^: a Stoner model [28 



a 



decimation method, a cluster RG |29[ ], and a Bethe-Peierls solution of the classical model p9fl. It 
appears that different treatments of the model have lead to different answers. For instance, there is 
a well defined transition in the Stoner criterion and the Bethe-Peierls treatments p8| , |29[ . On the 
contrary, there is no transition in the decimation method where the Hamiltonian Eqs. [l], § is mapped 



onto a percolation model [28]. A possible approach to this problem would be to generalize the work 
in Ref. p5| : instead of considering the model Eqs. p], |2] with infinite range exponential interactions, 
it is possible to approximate the problem by considering the Voronoi lattice model with exponential 
interactions, in which each lattice site has a finite number of neighbors. This model is well suited 
for carrying out numerical simulations, and avoids the difficulty that the initial model Eqs. |], || has 
infinite range interactions. 

Here, we would like to follow a different route and replace the original Hamiltonian Eqs. [l], |2] by 
a simplified one. This is done by noticing that the essential feature of the Hamiltonian Eqs. [|, |2] is 
that the exchange between two spins at a random position is distributed according to 

where we used the spacing distribution P(l) = xexp(— xl) and the exchange J(l) = Aexp(— //£). 
We replace the original model Eqs. [l], ^ by another model in which the spins are on the sites of a 
regular square lattice, and have a random nearest neighbor exchange in the distribution (|8|). Because 
the square lattice model with the exchange distribution Eq. || and the original model Eqs. [I], § are 
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E E 

Figure 6: The notations used in the renormalization. 

controlled by the same type of disorder, it is natural to conjecture that the two models have an 
identical physics. One way to study the square lattice model with the exchange distribution Eq. |8| 
would be to perform large scale Monte Carlo simulations. There is however a more direct way to 
handle the model, which consists in replacing the square lattice by a recursive hierarchical lattice 
(see Fig. ||), where the Migdal Kadanoff RG equations can be obtained in an exact form. We will 
show that a non trivial physics is going on in the hierarchical lattice model, which is an indication 
there is also a non trivial physics in the square lattice model with the exchange distribution Eq. 0. 



3.2 Renormalization group equations 

Let us derive the RG equations of the Ising model with the exchange distribution Eq. |8] on a hierar- 
chical lattice. The partition function associated to the exchange configuration on Fig. 0-(a) reads 



Z(E, £') = ^ exp /3 [JiEct + J 2 Za' + J 3 SV + J 4 EV] 



(9) 



and we impose that Eq. ^ be identical to the partition function associated to the exchange configu- 
ration on Fig. |||-(b), up to a proportionality factor: 



Using the relation 



we find J = Ji_3 + J2-4, where 



Z(E,E') = A/" exp (/?«/££'). 
W |Z(+,-) 



Jl - 3 = ^ ln 



cosh( / 3(J 1 + J 3 )) 
cosh(/?(Ji- J 3 )) 



J: 



2-4 



2(5 



cosh(/3(J 2 + J 4 )) 
cosh(/3(J 2 - J 4 )) 



(10) 



Now that we have determined the renormalization of the exchanges, we iterate the exchange distri- 
bution. Noting P{J\) ... P(Ji) the distribution of the exchanges J\ ... J 4 , and P(J) the distribution 
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x=0.02 




Number of RG iterations 

Figure 7: RG flow of the classical model relevant for the antiferromagnetic transition in CuGeOs. The 
parameters are A = 44.7 K and £ = 10, and we used a cut-off N = 300. The doping concentration 
is x = 0.02. The solid line represents the evolution of the average bond length ((/)). The doted 
lines represent the evolution of the width of the bond length distribution y/ ({(I — {(l))) 2 ))- The RG 
trajectories with T = 5 (+), T = 4 (x), T = 3 (*) and T = 2 (□) flow into the paramagnetic phase 
(the bond length renormalizes to infinity). The trajectories with T = 1 (■), T = 0.26 (•), T = 0.06 
(a) flow into the ferromagnetic phase (the bond length renormalizes to zero). 



of the renormalized exchange J, we have 



P{J) = / dJidJ 2 dJ3dJ4P(Ji)P(J2)P(J 3 )P(J4)6 



J — Jl_ 



1-3 



J: 



2-4 



(11) 



It will be useful to change variables to the bond lengths I = £ln (A/ J) and iterate the bond length 
distribution p(l) instead of the exchange distribution P{J). 

To calculate numerically the iteration of the exchange distribution Eq. [ll], we use a discrete bond 
length / = 1, ...,N, and we introduce an upper cut-off for the bond length. This is valid if we can 
check that the RG flow does not depend on N. 

3.3 Analysis of the RG flow 

The RG flow of the model is shown on Fig. [?]. We have shown on this figure the evolution of the 
average bond length ((I)) and the width of the bond length distribution \J (((I — ((l))) 2 ))- At low 
temperature, the average bond length renormalizes to zero (ordered phase) while at high temperature 
it renormalizes to infinity (paramagnetic phase). Therefore, there is a well-defined transition in the 
model. We have checked that the transition temperature is independent on the cut-off N used to 
iterate the bond length distribution. The existence of a thermodynamic transition could have been 



10 



anticipated on the basis of the simplest possible approximation of the RG flow (see Appendix |X|) . 
What is less obvious is that, after a transient in the first RG iterations, the width of the exchange 
distribution becomes much smaller than the average exchange: in spite of a broad initial exchange 
distribution in which ((/)} = yj (((I — ({I))) 2 }), the system renormalizes to an almost disorder- free 
exchange distribution in which y (((I — {(l))) 2 }) <C ((0)- Therefore, at large scale, the spin system 
looks ordered while inhomogeneities are visible only at small scale. This type of behavior may explain 
why there is a pronounced maximum in the temperature dependence of the susceptibility even at 
very low doping |lj 



4 Nature of the quantum disordered region 

Now, we would like to investigate the behavior of the model in the quantum disordered region, and 
compare it to experimental data on the Haldane gap compound Y2BaNiOs. More specifically, we 
would like to determine whether the behavior of the model is compatible with the susceptibility 



experiments by Payen et al. [18], who reported that the susceptibility of Y2BaNiOs has a power-law 
temperature dependence. In our opinion, it is an important question to determine which ingredients 
should be incorporated in the theoretical model to describe the existing experiments. A first strategy, 



followed by Batista et al. [46] is to look for the "most realistic possible" model. The first step in this 



approach is to consider that the relevant Hamiltonian for Y2BaNiOs takes the form 

H = V { JS t .S l+1 + D(S* f + E [(Sff - (S* ) 2 ] } • (12) 



The anisotropy parameters in Eq. 12 have been determined by fitting inelastic neutron scattering 



experiments jl7|, Next, a density matrix renormalization group (DMRG) method has been used 



to treat the Hamiltonian Eq. 12 in the presence of magnetic impurities. The authors of Ref. J46| are 
then able to reproduce specific heat experiments, and also arrive to an agreement with susceptibility 
experiments [ p| . The main objection that one might be tempted to formulate is that the Hamiltonian 
of the spin-1 chain relevant for Y2BaNiOs contains already three adjustable parameters, and that two 
more additional gyromagnetic factors have been added to describe the susceptibility experiments [18|. 



Here, it is proposed that the power law temperature dependence of the susceptibility is a generic 
feature of the quantum disordered region of the phase diagram, that can be explained in a model 
with only a minimal number of ingredients. 

4.1 Two-spin model 

As we already explained in the Introduction, in the low energy model relevant to describe the doped 
Haldane gap compound, each paramagnetic impurity generates a unit of two "edge" spin-1/2 moments 
(see Fig. |l|). The two spin-1/2 moments in the same unit are coupled by a ferromagnetic exchange 
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Figure 8: Log-log plot of the temperature dependence of T\{T) for the two-spin model and the cluster 
RG of the full chain. We used A = 100 K, £ = 6. For the two-spin model, we used x = 0.01 (□), 
x = 0.03 (©), x = 0.05 (A), x = 0.07 (V). The same symbols filled in black have been used for the 
cluster RG calculation. The product T\{T) has been fitted to a power-law dependence T\{T) ~ T a , 
with a = 0.017 {x = 0.01), a = 0.037 (x = 0.03), a = 0.04 (x = 0.05), a = 0.04 (x = 0.07). 
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Figure 9: Magnetization of the two-spin quantum model, with the parameters relevant for Y2BaNi05: 
A = 100 K, £ = 6. 
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J 2 the magnitude of which is of order of interchain interactions [44]. As shown on Fig. |], there is a 
staggered exchange (— )' exp (— //£) coupling two edge moments at a distance I. 

Let us first consider the simplest model in which the impurities are assumed to cut the chain into 
finite segments: J2 = 0. This is a valid model if the temperature is larger than J2, or equivalently, 
than the strength of interchain interactions. Consider two edge spins at distance I, coupled with a 
Heisenberg Hamiltonian H = J(7)Si.S2 — h(S\ + S^), and the exchange J(l) = A(— )' exp (— Z/£), 
with V(l+) = xexp(— x/+) the distribution of even length segments and V(l-) = xexp(— xl-) the 
distribution of odd length segments. It is easy to calculate the average magnetization in a magnetic 
field h: 

((M(h,T))) 



1 f A 

-x / P(J)[M+(h,T,J) + M-(h,T,J)]dJ, 

2 Jo 



with PdJl) = (x£/A)(| J|/A) 1+x ^ the exchange distribution, and 



M+(h,T, J) 



,(5h 



-f3h 



M_(/i,T, J) 



-Ph 



(13) 



(14) 



e ph + e -ph + i + e pj ' — v> - ' "i e ph + e -ph + 1 + e -pj 
the magnetization of the spins coupled by an antiferromagnetic or ferromagnetic exchange J at a 



finite temperature T = 1/(3. To calculate the magnetization Eqs. 13, 14, it is convenient to integrate 
by parts: 

((M(h,T))) 



, T 

xA + x — 
A 



x£. rA/T 



u^f(u,h/T), 



with 



A 



a ph 



-Ph 



Ph 



-Ph 



+ 



e Ph + e -ph + l + e A/T e ph + e -ph + l + g-A/T 



/(«, 



l__d_ 

2d« 



-/3h 



-/3h 



+ 



e /3/i + e -Ph + 1 + e « e /3ft + e -/3/i + 1 + e -u 



(15) 

(16) 
(17) 



The resulting susceptibility is shown on Fig. ||, where it is visible that the product T\{T) ~ T a 
has a power-law behavior at low temperature, like what has been observed in experiments on 
Y2BaNi05 [18|. The magnetization in an applied magnetic field is shown on Fig. ^, and is close 
to the experimental observation This shows that this model with only two spins contains much 
of the physics as far as the temperature is above the interchain coupling Jj_. 



4.2 Cluster RG 

4.2.1 Quantum disordered phase I 

Now, let us consider the cluster RG of the model with a finite J2, and first consider the behavior 
of the model when the temperature is larger than J2. I refer the reader to Refs. [2^, [2^| for an 
explanation of the method, and just present here the results. First, if the temperature is above J2, 
we find that the exponent of the power-law Curie constant T\{T) ~ T a is identical to the one of 
the two-spin model (see Fig IsT) . This is not a surprise because above J2 the edge moments in the 
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Figure 10: Temperature dependence of the Curie constant, with the parameters relevant for 
Y 2 BaNi0 5 : A = 100 K, f = 6, and x = 0.03, x = 0.05, and x = 0.07. We renormalized a 
chain with 10000 sites and averaged over 100 realizations of disorder. We used J 2 = 0.3 K. 



same unit remain uncoupled and the cluster RG contains the same physics as the two-spin model. 
Moreover, the two-spin model is an exact treatment while the cluster RG is approximate. The 
agreement between the two methods shows the validity of the cluster RG method. 



4.2.2 Quantum disordered phase II 

Now, I consider the cluster RG of a single chain at a temperature smaller than J 2 



J ± . The 



temperature dependence of the Curie constant is shown on Fig. 10. It is visible that the Curie 
constant increases strongly with decreasing the temperature below J 2 . This is because when T ~ J 2 , 
the survival spin- 1/2 moments in the same unit are frozen into spin-1 moments. This phenomenon 
is not quantum in nature because it occurs also in the classical disordered model, which is analyzed 
in details below. 

Once the freezing into spin-1 units has been done at the energy scale J 2 , the resulting effective 
model is again the one of spin-1 objects with random exchanges, being ferromagnetic or antiferro- 



magnetic. This results in the low temperature quantum disordered phase II on Fig. 10. 

In the presence of a finite interchain coupling, one still expects the appearance of two types of 
quantum disordered regions. It is also expected that the quantum disordered phase II on Fig. 10 is 
a 3D random singlet state, with singlet formation between spins belonging to different chains. The 
cross-over to a 3D regime is however not expected to change the shape of the temperature dependence 
of the susceptibility (see Fig. [K|). 
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Figure 11: Temperature dependence of the Curie constant T\{T) in the classical analog of the doped 
spin-1 chain. The parameters are £ = 6, A = 100 K, J2 = 0.3 K, relevant for Y2BaNiOs. The Curie 
constant crosses over from 2x at high temperature to 4x at low temperature. 

4.3 Classical disordered model 

Now, let us determine to what extend the physics of the classical magnet resembles the physics of 
the quantum magnet. The Ising chain Hamiltonian reads 



with the Ising variables Gi corresponding to the consecutive edge moments, and Jj^+i determined 
according to the rules on Fig. |l| The even bonds correspond to a ferromagnetic exchange — J2 and 
the odd bonds correspond to an exchange J (I) = A(— ) exp [— with the spacing I drawn in the 
Poisson distribution P(l) = 2 exp (— xl). 

As shown in Appendix |B|, the Curie constant crosses over from 2x at high temperature to Ax 
at low temperature, when the temperature decreases below J2. Therefore, in the classical model, 
the product T\{T) increases monotonically with a decreasing temperature, while the opposite is 
observed experimentally in Y2BaNiOs [jl8|]. This is not unexpected because the classical model 
cannot describe the gapless Haldane phases. Nevertheless, the increase in the susceptibility below J2 
is very reminiscent of the behavior of the quantum model in the same temperature range. This is 
because the freezing of the edge moments of the same unit in a ferromagnetic alignment occurs both 
in the classical and quantum models. 

The scaling function of the magnetization can be calculated easily by iterating numerically the 
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Figure 12: Scaling of the magnetization M(h, T) = T^^G{x, h/T) in the classical analog of the doped 
spin-1 chain. We have shown G(x, h/T) as a function of h/T. Each curve corresponds to values of T 
varying from 0.9 to 3 and h varying from to 1. We used the parameters £ = 6, A = 100 K relevant 
for Y 2 BaNi0 5 . We used 7 ~ 0.05. 



magnetization distribution Eq. 20 and using the relation 

V{M, h = 0) e P hM 



V{M, h) 



J £ M ,P(M>,h = 0)eP hM ' 



to obtain the magnetization distribution in a finite magnetic field. The magnetization takes the form 
M(h, T) = T _7 G(x, h/T) in a given temperature range where the susceptibility can be approximated 
by x(T) ~ j 1 -!-7 # The scaling function is shown on Fig. 12, where it is visible that it has qualitatively 
the correct behavior in spite of the exponent 7 having the wrong sign compared to experiments. 
Therefore, the shape of the scaling function of the magnetization does not appear to be a crucial test 
to the model. 



5 Conclusions 

To conclude, the present work was intented to describe doping a spin-Peierls and a Haldane gap state 
in a unified framework. We have first shown how the relevant energy scales in the problem could be 
calculated from the analysis of a two-spin model, which allowed to discuss the phases of the model as a 
function of the doping concentration and interchain interactions. In the relevant temperature window, 
there are two distinct regions depending on the doping concentration and interchain interactions: (i) 
an antiferromagnetic region; (ii) a quantum disordered region, remarkably, this type of phase diagram 
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compares well with the known behavior of the spin-Peierls compound CuGeOs and the two Nickel 
oxides PbNi 2 V 2 8 and Y 2 BaNi0 5 . 

Next, we used our approach to investigate in more details the two possible phases of the model. 
We have shown that the physics in the antiferromagnetic region of the phase diagram is classical in 
nature and therefore we were lead to study the corresponding Ising model. We have replaced the 
original Hamiltonian Eqs. [l], |2] by another Hamiltonian having the same features, and presented the 
solution of the latter Hamiltonian on a hierarchical lattice structure. Interestingly, we find that the 
renormalized problem is non disordered. This is in agreement with the presence of a well-defined 



cusp associated to antiferromagnetism in the susceptibility [14]. 

In the "quantum disordered" region of the phase diagram, the physics is strongly controlled by 
quantum fluctuations. Already in a model with two spins only does the susceptibility have a power- 
law temperature dependence, very similar to the experimental observation [i~8[| . We have suggested 
that there is no need to introduce many coupling constants to reach an agreement between the model 
and experiments. Going beyond the level of a two-spin model, we have found the existence of another 
quantum disordered phase at low temperature. We have analyzed the behavior of the classical model 
and shown it contains already a physics relevant to the quantum model. 

A Projection of the RG flow on a trial exchange distribution 

We would like to present the simplest possible approximation of the RG equations obtained in 



section 3.2 in which we project the RG flow onto the single parameter distribution Pn(]^) — — L n ^. 
The initial distribution is obtained via the relation J pq(1)cII = J xexp (— xl)dl, withpo(0 = S(l — Lq). 
This leads to 

L^On(l±f). (18) 



Next, we start from the distribution p n (l), make one RG transformation Eq. 11, and impose that 
the iterated distribution has the same first moment as p n +i(L), from what we can determine the 
parameter L n +%: 

/3A 

Ln+1 = € ln In cosh [2/3A exp (-£„/£)] ' (19) 
Since we want to discuss the stability of the paramagnetic phase, we consider Eq. |l9| in the limit 
of a large bond length: L n+ i ~ — £ln(2/3A) + 2L n . In this limit, we find L n = £ln(2/3A) + 
(Lq — £ ln (2/3 A)) 2 n . Using Eq. we obtain a phase transition at the temperature 

T = 2Axg 
c l + xf 

The transition temperature is in agreement with what has been found in previous works with different 



methods (the Stoner criterion [28| and the Bethe-Peierls method [|29f| ). As we show in the body of 
the article, the renormalized exchange distribution can be well approximated by the distribution 
Pn(l) = $(1 ~ Ln) in the sense that the problem renormalizes to a non disordered one. 
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B Solution of the classical analog of the doped spin-1 chain 



We consider a finite chain with N edge moments in which the end spin at site N is frozen in 
the direction +, and note V%{M) the corresponding magnetization distribution. We note Xi = 
exp (/?Jj)/[exp (PJi) + exp(— /3Jj)] the probability to find the spins Oi and cjj + i in an antiparallel 
alignment. We have 

V+ +1 (M) = (1 - x N )V+{M - 1) + x N Vjf{M - 1). (20) 

Using the relation V^{M) = Vjf(-M), we get 

(M}+ +1 = (l-2x N )(M)++l (21) 
(M 2 ) N+1 = (M 2 ) N+1 + 2(l-2x N )(M)++l. (22) 

These relations can be solved analytically. For this purpose, let us separate the even bonds cou- 
pled by the ferromagnetic exchange — J2 and note xf = e~@ j2 / [e^ j2 + e~^ j2 ]. The odd bonds are 
ferromagnetic or antiferromagnetic and we are lead to define 

e -/3j(o 

y = /-j V ^ ePJ(i) +e -pj(i)' 
The magnetization of a chain with an even number of sites N = 2p is found to be 

<M>+ = 2i lZx F) [1 " XP \ > (23) 

with X = (1 — 2xf)(1 — 2y). The correlation length is £t = —2/lnX. With an odd number of sites 
N = 2p + 1, the magnetization is 

The expressions of the first moment Eqs. |23], are next used to solve for the second moment: 

{M * hr ~ 2 P {i + ^-W-'ri + a(i-2^Ki-») 

to leading order in the chain length 2p. At high temperature, one has y ~ 1/2, xjr ~ 1/2 and X ~ 0, 
and therefore a susceptibility scaling like x ~ 2x/T. At low temperature, one has ~ 0, y ~ 1/2 
and X ~ 0, and a susceptibility scaling like x ~ 4x/T. Therefore, the Curie constant crosses over 
from the high temperature value 2x to 4x at low temperature, because the spin-1/2 moments in the 



same unit are frozen ferromagnetically at a temperature T ~ J2 (see Fig. 11 
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